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Abstract 


Desorption  of  an  adatom  from  a  solid  surface  through  vibrational 
excitation  by  an  IR  laser  is  studied  using  a  generalized  master  equation 
approach  and  the  isomnesic  approximation.  A  sharp  transition  from  no 
desorption  to  almost  instantaneous  desorption  is  seen  as  the  bound-to- 


continuum  population  transfer  rete  is  veried.  The  critical  rate  increases 
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I.  Introduction 

In  the  field  of  laser -stimulated  surface  processes  (LSSP),  desorption  Is  the 

most  commonly  studied  phenomenon.  As  demonstrated  In  a  number  of 

experiments,^  ^  the  desorption  cross  section  can  be  measured  unequivocally, 

albeit  with  some  difficulty.  One  aspect  of  these  experiments  which  Is 

especially  significant  Is  the  strong  dependence  of  the  desorption  signal  on  the 

laser  frequency.  The  control  that  this  selectivity  gives  us  In  directing  the 

course  of  chemical  reactions  on  a  solid  surface  holds  the  promise  of 

establishing  a  totally  new  kind  of  chemistry.  However,  not  only  Is  there  a 

shortage  of  well-characterized  laser-desorption  experiments,  but  also  a  clear 

understanding  of  the  detailed  mechanisms  responsible  for  the  selectivity  has  been 

lacking.  The  reader  Is  referred  to  recent  reviews  for  details.^  ^ 

A-7 

The  aim  of  our  early  work  on  LSSP  was  to  delineate  the  coarse  features  of 
the  phenomena  Involved.  Using  models  In  which  the  adspecies  (a  molecule) 
vibrated  In  the  attractive  region  of  a  surface  potential,  the  degree  of 
selectivity  was  shown  to  depend  on  the  anhara»nlclty  of  the  potential  and  the 
mismatch  between  the  laser  frequency  and  the  fundamental  vibrational 
frequency  of  the  adspecies  In  the  adsorption  potential.  In  addition,  the 
existence  of  energy  feedback  was  made  plausible.^  Thus,  assuming  that  the  laser 
photons  were  Absorbed  directly  by  a  so-called  internal  active  mode  A  of  the 
molecule,  idilch,  in  turn,  was  coupled  to  sink  modes  B  (other  modes  of  the 
molecule  end  phonon  modes  of  the  solid),  energy  was  shown  to  return  from  B  to  A 
cn  s  short  time  scale  for  specific  ranges  of  laser  power  and  system  parameters. 
While  this  cculd  be  rstlonsllsed  In  terms  of  s  slow  energy  dissipation  rate 
Inside  the  solid,  leading  to  local  buildup  In  the  vicinity  of  the  adsorbed 
molecule.  It  was  not  possible  to  distinguish  specific  mechanisms  within  that 
phencsMnolcgloal  framework. 


3 


More  mechanistic  approaches  have  subsequently  been  adopted,  generally 

g 

based  on  Morse  potentials  for  the  adspecles.  Mechanisms  for  desorption  based 
on  pumping  by  Internal  mode  excitation  of  the  adsorbed  molecule  and  tunnelling 
between  degenerate  bound  and  continuum  levels  of  the  adspecles^^  have  explained 
some  features  of  the  experimental  results.  In  an  attempt  to  Isolate  the 
conditions  for,  and  onset  of,  selectivity,  we  have  presented  a  time-dependent 
approach.  Phenomenological  time-independent  coupling  constants  are  replaced  by 
time-dependent  memory  kernels,  and  time  averages  are  supplanted  by  temporally 
varying  distribution  functions.  So  far,  energy  transfer  has  been  studied  In  the 
absence  of  a  desorption  channel,  but  selective  absorption  has  been  confirmed  for 
vanishing  mismatch  In  these  first -principles  quantum  mechanical 

treatments,  Is  represented  by  a  transition  frequency  (where  S  and  S'  are 
vibrational  quantum  nusbers  of  stationary  states  of  the  adsorption  potential) 
closest  to  Furthermore,  the  time -dependent  profile  of  the  average 
vibrational  energy  of  the  adsorptive  bond  (the  adbond)  Miows  an  initial  strong 
increase  followed  by  a  pronounced  minimum  and,  subsequently,  another  increase. 
The  decrease  leading  to  the  mlninm  Is  due  to  energy  transfer  to  the  bulk 
phonons;  and  the  subsequent  increase  Is  due  to  feedback  from  the  phonons  to  the 
adbond.  In  this  paper,  desorption  Is  considered  in  the  form  of  a  transition 
from  an  excited  bound  state  of  the  adsorption  potential  to  the  continuum.  The 
resulting  formalism  is  found  to  depart  insignif ieantly  from  that  of  our  previous 
Mork^*’*^  in  terms  of  interpretive  power.  Vithout  recourse  to  massive  oomputer 
ealoulations,  it  is  possible  to  obtain  time -dependent  probibillty  and  energy 
profiles  within  the  isomneslo  approx  last  Ion.  Details  of  the  desorption 
fornslisa  are  presented  In  Section  III,  which  follows  a  brief  aumary  of  the 
time -dependent  theory  of  vibrational  energy  transfer  between  an  adatom,  a  solid 
and  an  IR  laser  in  Section  IX.  Results  are  presented  and  dlacuseed  in 


Section  IV. 


111.  De sorption 

Desorption  corresponds  to  transfer  of  population  from  the  manifold  of  bound 
levels  of  the  adsorption  potential  to  the  continuum.  He  Introduce  a 
phenomenological  constant  transition  rate  R  from  the  highest  bound  level  |o> 
attainable  through  the  combination  of  resonant  laser  and  phonon  excitations,  to 
continuum  states  |k>.  In  our  case,  |o>  -  |3>.  As  a  result  we  must  Include  an 
extra  term  -RP,(t)  In  the  GME  for  Mt),  so  that 

0  o 

-  I  <*t'CKg3,(t  -  f)  ^3) 


and 


•  I  '  RP.r(t) 


S'fa 


S' o' 


.(R) 


In  the  Isomneslc  approximation,  we  replace  Eq. (2)  by 

*  **SS*  * 

where  Q33,  and  are  constants  In  time.  He  can  write,  for  all  S, 

•  I  IKsS-  ‘0 


(5) 


S'fS 


'S'S  '0 


®s«»  "  ^®S'S  * 


«‘s,J  l>s(t)|  . 


(6) 


Differentiating  and  taking  the  Laplace  transform,  we  get 
8^^3(8)  -  sPs(O)  -  ?s(0)  -  I 


S'-S 


'SS'^S' 


S' 


(7) 


all  S 


•SS''"'S 


%<here  we  have  defined  "diagonal"  elements  of  the  matrices  k  and  0  as 


6 


®SS  "  "  “S'S  "  * 

and(Pg(s)  is  the  Laplace  transfora  of  P^Ct). 

For  the  previously  studied  case  of  three  bound  levels,  we  get 


-A  B 


A  -B 


where. 


0  -R 


A  -  flo«  . 

B  -  Qq{1  ♦  a)  , 
o  -  C'/2x  . 

X  .  A-p/2kgT  , 
C*  .  e"P*  , 


p  being  a  fitting  parameter  Involved  in  the  integration  over  the  phonon  density 
fi 

of  states,  kg  the  Boltzmann  constant,  and  T  the  temperature.  The  local  rate 
due  to  the  phonon-adbond  interaction  is  obtained  from  Qgg'*' 

«SS.  ■  «o 

Inversion  of  the  Laplace  transform  equation  Involves  the  solution  of  the 
polynomial  equation 


where 


q(s)  ■  det  M  -  0  , 


k  ♦  sO  -  s^A  , 


A  being  the  unit  matrix.  Thus  we  need  the  roots  of  the  polynomial  equation 


•(s  ♦  As) 


-(s  ♦  Bs  ♦  kg) 


-(s‘  ♦  Rs  ♦  kg) 
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where  •  kg^,  for  the  pair  of  levels  resonant  with  the  laser.  Equation  (20) 

is  a  quartic  for  R  ^  0  and  a  cubic  for  R  -  0.  The  roots  can,  in  principle,  be 

obtained  in  closed  form,  but,  because  of  their  complicated  nature,  are  of  little 

practical  use  in  our  analysis.  Ve  have  instead  adopted  a  numerical  approach 

whose  details  we  shall  omit.  The  resulting  probability  profiles  are  used  to 

define  the  desorption  probability  and  average  adbond  energy  as  follows: 

Pp(t)  -  P^(0)  -  I  Pg(t)  .  (21) 

s 

£(t)  -  |es  Ps  .  (22) 

where  is  the  energy  eigenvalue  corresponding  to  IS^,  end 

Pt(0)  -  I  Ps<0)  •  (23) 

The  value  of  P^(0)  is  slightly  smaller  than  unity  because  the  set  of  levels  {s} 
used  in  the  GHE  does  not  include  all  bound  states  of  the  effective  adbond 
potential. 

The  specific  nature  of  the  parameters  A,  B  and  is  of  interest  from  an 
analytic  point  of  view.  Their  form  can  be  derived,  frMD  exact  expressions  for 

g 

the  memory  kernels  kgg,  presented  in  earlier  work,  by  using  the  lA,  Eq.  (5).  In 
light  of  Eqs.  (12)'(16),  the  parameters  A,  B  and  Eq  involve  and  kgg,,  which 
are  given  by 

“o  "  %  ®p  *SS'  •“'SS'I  " 

and 

*'ss'  “  V‘ss'^^ 


» 


(25) 
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whers 


Ap  -  . 

Bp  -  6(p<0>k3T)2/M,Hu^  . 

"SS'  '  “s‘s  ®ss* 

w* 

1^?  -  <S|.  |S'>  , 

-  (I.J/cjCll^A<“>*)  . 


*SS'  " 


(26) 

(27) 

(28) 

(29) 

(30) 

(31) 


Thtt  quantities  in  Eqs.  (26)>(31)  involve  the  following: 

t  8^^^:  well  depth  and  steepness  paraaeters  for  an  effective  Morse 

potential  between  the  surface  atom  and  the  adatom 

N^:  atomic  mass  of  the  solid 

take  on  the  values  1  and  2 
~Z 

u  :  mean  square  displacement  of  the  atoms  of  the  solid  from  equilibrium. 

s:  position  normal  to  the  solid  surface  measured  with  respect  to  the  equilibrium 
position  of  the  outer  most  lattice  atom. 

Co) 

Sq  !  location  of  the  minimum  of  the  effective  Morse  potential  between  the 
surface  atom  and  the  adatom 
It  laser  intensity 

Oq:  charge  difference  between  the  adspecles  and  the  surface 
Cq!  permittivity  of  free  space. 
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The  strengths  of  the  adbond-solld  interaction  and  the  laser-adbond 
interaction  are  represented  by  and  k^.  respectively.  It  has  been  demonstrated 
that  the  two  must  be  comparable  in  order  for  the  system  to  manifest  significant 
synergistic,  selective  effects. 

IV.  Results  and  Diacusalon 

With  the  inclusion  of  a  unidirectional  desorption  channel,  we  expect  the 
total  population  of  bound  states  to  be  depleted  monotonlcally,  1. e. ,  P^Ct)  must 
grow  monotonlcally  to  a  maximum  possible  value  of  P^(0).  Similarly,  since 
desorption  is  due  to  transfer  of  population  from  the  high-lying  level  |o>  to  the 
continuum,  the  average  adbond  energy £(t)  might  be  expected  to  increase  with 
time.  However,  because  of  the  competition  between  laser  pumping,  which 


populates  the  level  |o>,  and  desorption,  which  depletes  it,  the  final  behavior 


of  £(t)  is  not  as  simple  as  that  of  Pjj(t).  For  vanishing  R,  £(t)  should  have 

A _ ... _ _ ®  10  _ _ 1.  _ _ ^  ^ 


the  familiar  nondesorptlve  profile. 


For  values  of  R  comparable  to  the 


effective  laser  pumping  rate  (characterized  by,  but  not  equivalent  to,  k^)  the 
behavior  of  £(t)  cannot  be  predicted  on  intuitive  grounds.  For  the  other 
extreme,  namely  for  laser  pumping  rates  much  smaller  than  the  rate  R, 

the  average  adbond  energy  is  expected  to  rise  quickly  as  population  is  depleted 
uniformly  fhom  the  lower-lying  levels  and  does  not  return  to  them.  Eventually 
the  system  should  be  completely  desorbed  and  £(t}  should  become  zero.  Implying 
the  complete  absence  of  populated  bound  states. 

Our  results  for  the  prototype  system  of  references  11-13  bear  out  these 
expectations.  Figures  2  through  10  show  the  desorption  probabilities  Pp(t)  and 
the  associated  average  adbond  energies  £(t)  for  laser  powers  ranging  ft*om 
10  W/cm  to  10  W/cm  and  for  desorption  rates  spanning  five  orders  of 
magnitude.  Notice  that  the  time  scales  over  which  the  system  exhibits 
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significant  dynamical  effects  Increase  with  decreasing  laser  power,  as  expected. 
Significant  dynamical  changes  take  place  more  slowly  for  lower  laser  powers. 

Both  Pp(t)  andfi(t)  display  a  transitional  dynamical  state  of  the  system  for 
a  certain  narrow  range  of  desorption  rates.  For  rates  below  this,  the  system  is 
relatively  tranquil,  with  the  laser  causing  almost  no  change  in  the  adsorbed 
system.  The  desorption  probability  changes  imperceptibly  over  quite  a  large 
time  period,  and£(t)  oscillates  about  an  eventual  steady-state  value.  The 
behavior  of  the  system  for  R  larger  than  this  transitional  range  is  radically 
different  and  violent.  The  adbond  energy  £,(t)  and  the  desorption  probability 
Pjj(t)  increase  rapidly  and  monotonlcally  to  new  equilibrium  values,  the  rise 
time  decreasing  rapidly  with  increasing  R.  The  changes  in  behavior  of  the 
system  with  I  and  R  are  quantified  in  Fig.  11,  where  regions  of  oscillatory  (0), 
nonosc illatory  (N)  and  transitional  (T)  behavior  are  demarcated.  For  the  range 
of  R  studied,  namely  10  ^  ps  ^  i  R  S  1  ps  ' ,  both  the  cessation  of  slow 
oscillatory  behavior  and  the  onset  of  fast  sigmoid  behavior  takes  place  at 
higher  laser  powers  for  larger  transition  strengths  R.  Vfhlle  the  actual 
locations  of  the  boundaries  between  the  regions  0  and  T  and  the  regions  T  and  N 

(obtained  by  inspection  from  Figs.  2-10)  are  crude  estimates,  the  width  of 
region  T  is  uniform  and  the  edges  essentially  linear.  The  trends  depicted  in 
Fig.  11  are  important  from  the  point  of  view  of  selective  laser  photochemistry. 
Thus,  for  a  given  value  of  R,  laser  powers  in  the  region  0  will  cause  a  very 
slow  accumulation  of  vibrational  energy  in  the  adbond,  while  laser  powers  in 
regions  T  and  N  will  excite  the  adbond  very  quickly  and  cause  almost 
instantaneous  desorption. 

The  actual  experimental  situation  is  complicated  by  the  fact  that,  in 
general,  there  is  a  contribution  to  R  due  to  the  phonon  field  itself.  The 
strong  interaction  between  the  adbond  and  phonons  would  therefore  dominate  the 
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desorption  dynamics,  and  sampling  a  range  of  R  to  test  for  the  two  types  of 
behavior  would  be  difficult.  For  the  case  of  a  bound  state  |6>  embedded  In  the 
continuum,  however,  resonant  excitation  from  |o>  to  |6>  using  another  laser 
provides  a  practical  means  of  varying  R.  A  fairly  sharp  Increase  In  the 
desorption  cross  section  as  a  function  of  laser  intensity  would  be  the  signature 
of  the  transition  from  a  mostly  bound  to  a  mostly  desorbed  state.  The  whole 
process  is  strongly  bond-selective,  and  for  the  case  of  a  variety  of  adsorbed 
species,  different  pairs  of  lasers  should  cause  different  adspecies  to  desorb. 
Also,  for  low-power  operation  with  lasers  whose  frequencies  do  not  match  system 
transitions,  no  sudden  transition  would  be  observed. 

Purely  thermal  effects  can  be  expected  to  dominate  for  the  case  of 
signifcant  departure  from  the  two  generic  resonance  conditions  shorn  in  Fig.  1. 
These  conditions  are  designed  to  highlight  the  mechanism  responsible  for 
selective  effects  associated  with  laser-stimulated  desorption.  An  additional 
factor  that  exerts  a  strong  influence  on  the  dynamics  of  LSSP  in  general  is  the 
rate  of  dissipation  of  vibrational  energy  from  the  Debye  phonons  to  the  rest  of 
the  allowed  phonon  states.  The  resonance  condition  involving  the  Debye  frequency 
translates  to  an  energy  transfer  rate  between  the  adbond  and  the  solid  which  is 
much  larger  than  the  phonon  dissipation  rate,  a  factor  of  central  importance  for 
the  results  of  the  model  presented  here.  The  implied  'localization*  of 
vibrational  energy  in  a  thin  slice  of  the  phonon  density  of  states  is  the  k-space 
auialog  of  spatial  localization  of  energy,  but  with  additional  many-body  effects 
exhibited  as  oscillations  of  £(t).  Thus,  for  example,  if  the  phonon  energy 
dissipation  rate  is  comparable  to  or  greater  than  the  laser  pumping  rate,  the 
'feed*  level  |3>  for  desorption  will  not  be  sufficiently  populated  to  exhibit  the 
large  increase  in  average  adbond  energy  due  to  selective  resonant  excitation.  On 
the  other  hand,  if  the  phonon  energy  dissipation  rate  is  comparable  to  or  greater 
than  the  energy  transfer  rate  between  the  adbond  and  the  Debye  phonons,  the 
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feedback  effects  represented  by  the  oscillations  of  £(t)  would  not  be  evident  on 
the  time  scales  of  Figs.  2-7.  The  occupation  of  Debye  phonons  would  not  get 
sufficiently  large  for  the  back  transfer  of  vibrational  energy  from  the  solid  to 
the  adbond  to  be  significant.  However,  within  the  framework  of  the  assusied 
resonance  conditions,  the  many-body  effects  lead  to  the  very  interesting 
phenomena  described  earlier.  Experimentally,  the  sharpness  of  the  resonances 
will,  of  course,  determine  if  the  selective  effects  are  observable  directly,  or 
if  they  will  be  obscured  by  nonresonant  thermal  effects.  It  would  be  interesting 
to  investigate  the  possibility  of  separating  the  selective  and  the  nonselectlve 
effects  in  the  experimental  results.  Knowledge  of  the  behavior  of  P|^  or  £  as  a 
function  of  R  for  a  fixed  time  (perhaps  an  asymptotic  limit)  due  to  purely 
thermal  effects  is  necessary  for  such  a  separation. 

The  GME  for  the  desorptive  case.  Cq. (6).  suggests  a  similarity  between  the 
treatment  of  the  phonon-adbond  interaction,  l.e.,  a  delta  function  limit  for 
Kg|,^(t),  and  that  of  the  desorptive  term.  While  there  are  two  contlnuua 
present  -  the  phonon  distribution  on  the  one  hand  and  the  desorption  states  on 
the  other  -  it  is  not  clear  that  the  two  can  lead  to  similar  memory  effects.  If 
the  I o>-to-continuum  transition  is  primarily  to  a  very  small  slice,  it  may  be 
necessary  to  treat  this  channel  analogously  to  the  laser-adbond  interaction, 
namely  with  an  Isomnesic  kernel.  A  detailed  treatment  with  well-defined 
excitation  mechanisms  is  necessary  for  a  better  understanding  of  the  role  of  the 
desorption  channel. 
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